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Abstract 

Interspecific hybridization is known for triggering genetic and epigenetic 
changes, such as modifications on DNA methylation patterns and impact on 
phenotypic plasticity and ecological adaptation. Wild potatoes (Solanum, 
section Petota) are adapted to multiple habitats along the Andes, and natural 
hybridizations have proven to be a common feature among species of this 
group. Solanum x rechei, a recently formed hybrid that grows sympatrically 
with the parental species S. kurtzianum and S. microdontum, represents an ideal 
model for studying the ecologically and evolutionary importance of hybridiza- 
tion in generating of epigenetic variability. Genetic and epigenetic variability 
and their correlation with morphological variation were investigated in wild 
and ex situ conserved populations of these three wild potato species using 
amplified fragment length polymorphism (AFLP) and methylation-sensitive 
amplified polymorphism (MSAP) techniques. We observed that novel methyla- 
tion patterns doubled the number of novel genetic patterns in the hybrid and 
that the morphological variability measured on 30 characters had a higher cor- 
relation with the epigenetic than with the genetic variability. Statistical compari- 
son of methylation levels suggested that the interspecific hybridization induces 
genome demethylation in the hybrids. A Bayesian analysis of the genetic data 
reveled the hybrid nature of S. x rechei, with genotypes displaying high levels 
of admixture with the parental species, while the epigenetic information 
assigned S. x rechei to its own cluster with low admixture. These findings sug- 
gested that after the hybridization event, a novel epigenetic pattern was rapidly 
established, which might influence the phenotypic plasticity and adaptation of 
the hybrid to new environments. 



Introduction 

Hybridization and polyploidization are important evolu- 
tionary forces in plants, as a source of genetic and epige- 
netic variation. Interspecific hybridizations generate genetic 
alterations and modifications on DNA methylation 
(Salmon et al. 2005; Marfil et al. 2006). The latter is con- 
sidered an epigenetic mark, along with histone modifica- 
tions and regulation by small RNAs. Epigenetic variations 
generate new phenotypic variants that could be targeted by 
selection (Kalisz and Purugganan 2004; Rapp and Wendel 
2005). Also, in many cases, epigenetic marks can be trans- 
mitted across generations operating as genuine genetic vari- 



ation, and therefore, they could play a significant role in 
speciation and evolution (Jablonka and Raz 2009). 
Amplified fragment length polymorphism (AFLP) and 
methylation-sensitive amplification polymorphism (MSAP) 
analysis showed that hybridization more than polyploidiza- 
tion triggered epigenetic changes in the allopolyploid 
genomes of Spartina (Ainouche et al. 2004), Senecio 
(Hegarty et al. 2011), and Arabidopsis (Wang et al. 2006). 
In recent years, some studies assessing the importance of 
epigenetic variation in natural populations were published 
(Salmon et al. 2005; Marfil et al. 2009; Paun et al. 2010; 
Lira-Medeiros et al. 2010), but it is necessary to increase 
empirical evidence for natural epigenetic variation. 
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Figure 2. Distribution of Solarium kurtzianum and Solarium microdontum and sampled regions in this work. (A) Distribution of S. kurtzianum 
(red, Argentina) and S. microdontum (blue, Argentine and Bolivia). (B) Map of La Rioja province and detail of two regions where natural 
populations of S. kurtzianum, S. microdontum, and Solanum x rechei were collected (yellow). 



Wild potatoes (Solanum, section Petota) constitute an 
ideal model to analyze the evolutionary consequences of 
hybridization in the generation of epigenetic variation and 
its significance in the phenotypic plasticity and ecological 
adaptation: (i) Natural hybridization and hybrid specia- 
tion have been proposed to be prevalent in this group of 
plants (Hawkes 1990; Masuelli et al. 2009); (ii) wild 
potato species are adapted to multiple habitats along the 



Andes in America and represent a group with extensive 
phenotypic plasticity (Hawkes 1990); and (iii) both modes 
of reproduction in wild potatoes, sexual by seeds and asex- 
ual by stolons and tubers, would allow the study of trans- 
generational inheritance and mitotic stability of epialleles. 

Analyzing synthetic hybrids between Solanum tubero- 
sum and Solanum kurtzianum, we observed DNA methyl- 
ation changes induced by hybridization (Marfil et al. 
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2006). In the natural hybrid Solatium ruiz-lealii, a rela- 
tionship between methylation patterns and abnormal flo- 
ral morphology was observed (Marfil et al. 2009). In this 
study, we expect to extend the observations described 
above by analyzing natural populations of an interspecific 
hybrid and its parental species. 

Of 27 taxa of putative hybrid origin within section 
Petota (Spooner and van den Berg 1992), Solarium x rec- 
hei (Hawkes et Hjerting) is one of the most extensively 
documented (Fig. 1). The distribution of S. x rechei is 
restricted to the locality of Guanchin, at the southern end 
of the Sierra de Famatina in Argentina, province of La 
Rioja, department of Chilecito. In this locality, the paren- 
tal species S. kurtzianum (Bitter et Wittm.) and Solanum 
microdontum (Bitter) also grow sympatrically (Fig. 1). 
Solanum kurtzianum is a diploid species (2n = 2x = 24) 
that grows on four Argentinean provinces, along a tran- 
sect approximately 900 km long (Hijmans et al. 2002; 
Fig. 2). It is adapted to dry regions within the Monte 
phytogeographical province (Cabrera 1976) with an 
annual average temperature of 13°C and an annual pre- 
cipitation of 314 mm (Hijmans et al. 2002). The distribu- 
tion range of S. microdontum extends from northern 
Argentina to southern Bolivia over a distance of about 
1300 km (Hijmans et al. 2002; Fig. 2). It grows typically 
in the rain forest, corresponding to the Cabrera's (1976) 
phytogeographical province of the Yungas, although it 
extends further south than the limits of the rain forest in 
La Rioja province, where it grows in wet shady places 
with an annual average temperature of 14.5°C and an 
annual precipitation of 639 mm (Hijmans et al. 2002). 
Most of the populations of S. microdontum are diploid, 
although in the extreme southern end of its range, which 
corresponds to the locality of Guanchin, only triploid 
(2m = 3x = 36) populations have been recorded (Okada 
and Hawkes 1978; Okada 1981). Solanum x rechei has a 
very narrow distribution, and it seems to be better 
adapted to dry places with an annual average temperature 
of 12.8°C and an annual precipitation of 242 mm 
(Hijmans et al. 2002). 

Solanum x rechei is a recently formed hybrid showing 
instabilities such as plants with severe floral abnormalities 
and low pollen fertility. Its hybrid origin was first sug- 
gested by Hawkes and Hjerting (1969). Okada and 
Hawkes (1978) further tested its hybridity by field stud- 
ies, reciprocal synthetic reconstructions of the hybrid, 
analyses of morphological characters, chromosome 
counts, and pollen stainability. Clausen and Spooner 
(1998) confirmed the hybrid origin of S. x rechei per- 
forming RFLP analysis on accessions of the hybrids and 
the two parental species. The aim of this work was to 
study the basis of the novel variability originated by a 
recent interspecific hybridization event. For this purpose, 



N. Cara ef al. 

we analyzed the epigenetic, genetic, and morphologic var- 
iability of wild populations and germplasm bank main- 
tained accessions of three sympatric Solanum species, the 
natural hybrid S. x rechei, and its parental species S. ku- 
rtzianum and S. microdontum. We found that morpho- 
logic variability correlated better with epigenetic than 
with genetic variability. Additionally, a Bayesian approach 
showed that S. x rechei individuals established novel epi- 
genetic patterns, different from those of the parental spe- 
cies, whilst maintaining genetically admixed genomes. 
These findings suggest that epigenetic variation may play 
an important role in the stabilization and evolution of 
this hybrid. 

Materials and Methods 
Plant material 

The study was performed on two types of materials: (i) 
Plants from wild populations of S. kurtzianum, S. micr- 
odontum, and S. x rechei sampled as tubers in 2008 and 
2010 at the department of Chilecito, La Rioja; and (ii) 
accessions of these species provided by the Potato and 
Forages Germplasm Bank Instituto Nacional de Tecno- 
logfa Agropecuaria, Balcarce, Argentina (PFGB-INTA; 
Fig. 2 and Table 1). The study was carried out over two 
years. In the first year, tubers collected by the authors 
were sown in pots, and seeds provided by the PFGB were 
germinated in Petri dishes and then transplanted into 
pots. In the second year, all plant materials were grown 
from the tubers obtained in the previous season; thus, all 
genotypes were bred at the same time under uniform 
conditions in an insect-proof screen house. 

Morphological analysis 

Thirty qualitative and quantitative characters were mea- 
sured in three plants per accession once flowering had 
started (Table 2). Morphological measurements were per- 
formed on the seventh true leaf from the base of the plant 
and on the first inflorescence to appear. To test for 
hybridity based on morphological intermediacy, the 
"character count procedure" was performed as described 
by Wilson (1992). This analysis consists in tabulating 
intermediate and non-intermediate character states in the 
putative hybrid respect to its parents. To determine 
whether differences in characters were statistically signifi- 
cant, univariate analysis of variance (ANOVAs) were 
performed for each quantitative character and chi-square 
tests were employed for qualitative characters using 
STATGRAPHICS Centurion XVI software (StatPoint 
Technologies, Warrenton, VA). Only characters that dif- 
fered between the parents were taken into account, and a 
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Table 2. Morphological characters measured in Solarium kurtzianum (ktz), Solarium microdontum (mcd), and Solarium x rechei (rch) and test for 
intermediacy of S. x rechei. 

Mean (SDVMedian 1 ANOVA/chi-square test 2 

Morphological character ktz mcd rch ktz-mcd ktz-rch mcd-rch Intermediacy 
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Significance levels: ** = P < 0 .01; * = P < 0.05; ns = P > 0.05. 

1 For quantitative characters (1-21), mean and standard deviation (SD) were calculated. For qualitative characters (22-30), the median was 
obtained. 

2 Univariated ANOVA tests were performed for quantitative characters and chi-square tests for qualitative characters. 



one-sided sign test (Zar 1984) of intermediate versus 
non-intermediate characters was applied to judge whether 
the coalescence of intermediate character states was too 
improbable to represent divergence in the same characters 
and in the same direction. /Also a principal coordinates 
analysis (PCoA) was carried out, and a minimum span- 
ning tree (MST) was superimposed. As we had qualitative 
and quantitative variables, the Gower dissimilarity coeffi- 



cient (Gower 1971) was used to generate a distance 
matrix in the PCoA. The operational taxonomic unit 
(OTU) was the average of three plants. A cluster analysis 
was performed using the same distance matrix and UP- 
GMA (Unweighted Pair-Group Method with Arithmetic 
Averaging) as linkage method. R software 2.15.1 (R Core 
Team 2012) was used to perform the one-sided sign test 
in the character count procedure and to obtain the Gower 
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dissimilarity matrix. The NTSYS-pc 2.11 software (Rohlf 
1998) was employed for PCoA and cluster analysis. 

Amplified fragments length polymorphism 
(AFLP) analysis 

DNA was extracted from fresh leaves following Dellaporta 
et al. (1983) with an intermediate treatment with RNase. 
After quantification of DNA concentration by spectropho- 
tometry (GeneQuant RNA/DNA Calculator, Pharmacia 
Biotech), samples were diluted to a concentration of 
100 ng-/iL~\ 

To assess genetic variability in the evaluated material, 
an AFLP analysis was performed as described by Vos 
et al. (1995). £coRI and Msel were used as the rare and 
frequent cutter, respectively. The restriction reaction con- 
tained 2 units of EcoRI, 2 units of Msel, 1.25 /(L of buffer 
2 (NEB), 100 ng-/jL _1 BSA, and 250 ng of DNA in a final 
volume of 12.5 /iL and was incubated for 3 h at 37°C. 



Selective primer combinations (£coRI+3 and Hpall - Mspll+3) 
EcoRI + AGA/Hpall - Mspl + ATC 
EcoRI + ACA/Hpall - Mspl + AAT 
EcoRI + ACA/Hpall - Mspl + ATC 
EcoRI + AAC/Hpall - Mspl + ATC 
EcoRI + AAC/Hpall - Mspl + AAT 

Total: 
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The ligation reaction contained 5 pmol of EcoRI adaptors, 
50 pmol of Msel adaptors, 1.25 /iL T4 DNA ligase buffer 
(Promega), 0.75 units of T4 DNA ligase (Promega, Madi- 
son, WI), and 6.25 /iL of the digested products in a final 
volume of 12.5 /jL. Adaptors were ligated for 3 h at 20°C. 
Preselective amplification was accomplished using 
[£coRI + A] and [Msel + C] primers. Each PCR con- 
tained 40 ng of each primer, 400 /miol of dNTPs, 1 x 
Taq DNA polymerase buffer, 1 unit of Taq DNA poly- 
merase (Invitrogen), and 1 /(L of the ligation products in 
a final volume of 25 [iL. The amplification profile was as 
follows: 20 cycles at 94°C for 30 s, 56°C for 1 min, and 
72°C for 1 min. The pre-selective PCR products were 
diluted 1:3 with sterile ultra-pure H 2 0. Selective amplifi- 
cation was performed with fluorescence-labeled £coRI- 
selective primers and unlabeled Msel selective primers 
(Table 3). The reaction contained 0.125 /miol-L~ of 
[Msel + 3] primers, 0.125 /miol-L -1 of [£coRI + 3] 
labeled primers, 200 jUmol-L" 1 of dNTPs, lx Taq DNA 



Number of fragments 



137 
64 
197 
191 
73 
170 
832 



46 
46 
43 
40 

57 
232 
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Table 3. AFLP and MSAP adaptors and primers used in this study and fragments obtained. 



Sequences 



AFLP oligonucleotides 
Adaptors 
EcoRI 

Msel 

Preselective primers 
EcoRI + A 
Msel + C 

Selective primer combinations (EcoRI+3 and Msel+3) 
EcoRI + ACG/Msel + CTG 
EcoRI + ACG/Msel + CTC 
EcoRI + ACA/Msel + CTC 
EcoRI + AAG/Msel + CAC 
EcoRI + ACA/Msel + CTG 
EcoRI + AAG/Msel + CTG 

MSAP oligonucleotides 
Adaptors 
EcoRI 

Hpall - Mspl 

Preselective primers 
EcoRI + 0 
Hpall - Mspl + A 



S'-AATTGGTACGCAGTCTAC 
5'-CTCGTAGACTGCGTACC 
S'-TACTCAGGACTCAT 
S'-GACGATGAGTCCTGAG 

5'-AGACTGCGTACCAATTCA 
5'-G AC GATGAGTC CTG AGTAAC 



Total: 



S'-CTCGTAGACTGCGTACC 
S'-AATTGGTACGCAGTCTAC 
5--GACGATGAGTCTCGAT 
5'-CGATCGAGACTCATC 

5--GACTGCGTACCAATTC 
5'-ATGAGTCTC G ATC G G A 
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polymerase buffer (Invitrogen), 0.5 units of Taq DNA 
polymerase (Invitrogen), and 5 jiL of the diluted pre- 
selective product in a final volume of 20 /iL. The PCR 
protocol was as follows: 14 cycles at 94°C for 30 s, 65°C 
for 30 s (decreasing 0.7°C per cycle), and 72°C for 1 min, 
followed by 23 cycles at 94°C for 30 s, 56°C for 30 s, and 
72°C for 1 min. 

After selective amplification, 1 /iL of each sample was 
mixed with 8.5 /.tL Hi-Di formamide and 0.5 pL of the 
size standard Genescan 500HD ROX (Applied Biosystems, 
Foster City, CA). Samples were then denatured at 90°C 
for 3 min and were run on an Applied Biosystems 
Genetic Analyzer 3130. Semi- automated scoring was per- 
formed on the resulting electronic profiles using GeneM- 
apper version 3.7 (Applied Biosystems). 

Methylation-sensitive amplified 
polymorphism (MSAP) analysis 

Epigenetic variability was estimated by MSAP analysis, 
which characterize the methylation pattern at anonymous 
5'-CCGG sequences from random genomic DNA. The 
MSAP procedure was performed according to Xiong et al. 
(1999), with slight modifications. This is an adaptation of 
the original AFLP protocol, used to detect methylation 
polymorphism. It substitutes the frequent cutter enzyme 
Msel by Hpall and Msp\, while the rare cutter remains 
unchanged. Hpall and Mspl are isoschizomers which 
recognize the same tetranucleotide restriction site 
(5'-CCGG), but have different sensitivities to methylation 
at cytosines. Hpall is sensitive to full methylation (both 
strands methylated) of either cytosine but cleaves the 
hemimethylated external cytosine, whereas Mspl is sensi- 
tive only to methylation of the external cytosine. In this 
technique, the presence of the fragments in both profiles 
indicates a non-methylated site, while the absence of these 



fragments in both amplifications is ambiguous. It could 
be due to complete site methylation or to variation in the 
nucleotide sequence (Fig. 3). 

In the first restriction assay, 700 ng of genomic DNA 
was digested with 10 U of EcoRl (Promega) with 
100 ng-jUL -1 BSA and lx buffer H (Promega) in a final 
volume of 20 /(L. To ensure complete digestion, incuba- 
tion was performed at 37°C over night. The digest was 
split into two equal volumes and 10 U of Hpall 
(Promega) and 10 U of Mspl (New England Biolabs) were 
added to each half, along with 100 ng-/iL~' BSA and the 
appropriate buffers in a final volume of 20 /.iL for the sec- 
ond digestion. Ten /(L of the digested solution was ligated 
with 25 pmol of each EcoRl and Hpall/Mspl adaptors, 
using 0.75 units of T4 DNA ligase (Promega) in a final 
volume of 20 /(L. Adaptors were ligated for 3 h at 15°C. 
Pre-amplification was performed using 1 fiL of the liga- 
tion products and 4 pmol of the [EcoRl + 0] and [Hpall/ 
Mspl + A] primers in a final volume of 20 /(L containing 
lx PCR buffer, 0.1 mmol-L -1 dNTPs, and 1 U Taq poly- 
merase (Invitrogen). The pre-amplification products were 
diluted 1:3, and 1 fiL was used in the selective amplifica- 
tion reaction with the EcoRI- and Hpall/MspI-selective 
primers (Table 3) in a final volume of 20 ptL. The other 
components were the same as the pre-amplification reac- 
tions. In the pre-amplification as well as in the amplifica- 
tion, PCR programs were the same as in the AFLP 
protocol. The final amplification products were separated 
by electrophoresis on a 6% denaturing polyacrylamide 
gel, silver stained following Benbouza et al. (2006) and 
scanned for manual scoring. 

Data analysis 

In the AFLP analysis, a binary matrix was constructed 
scoring fragments as present (1) or absent (0). For the 
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Figure 3. Methylation-sensitive amplified 
polymorphism (MSAP) band patterning 
determined by HpaWMspi isoschizomers and 
correspondence with methylation status of 
CCGG sites. Black squares indicate methylated 
cytosines. 
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MSAP analysis, two different matrices were built: (i) a 
binary matrix with methylation-insensitive polymor- 
phisms, loci where the fragments were either present or 
absent in both EcoRl/Hpall and EcoRl/Mspl digests for all 
the genotypes; and (ii) a methylation-sensitive polymor- 
phisms matrix, with fragments that differed in the pres- 
ence/absence patterns between the EcoRl/Hpall and the 
EcoRllMspl digests in at least one genotype. In this case, 
patterns were codified from 0 to 3 (Fig. 3), and then, this 
codification was converted into a binary matrix for pres- 
ence (1) or absence (0) of the particular pattern. A gen- 
eral methylation level was calculated for each species by 
computing the percentage of methylated, hemimethylated, 
and unmethylated loci, determined by patterns 3, 2, and 
1, respectively. As absence of bands (pattern 0) could be 
due to methylation of external cytosines as well as varia- 
tions in the nucleotide sequence, they were not consid- 
ered as methylated loci. 

For both AFLP and MSAP procedures, a subset of 23 
samples were replicated starting from the same DNA 
extractions, to calculate the error rate of the techniques. 
Loci with different patterns on replicates were excluded 
from analysis as possible methodological artifacts. 

The total number of fragments (loci) for AFLP and 
MSAP was computed for each of the three species. Com- 
mon and species-specific fragments were identified, and 
the mean number of fragments per individual was also 
calculated. Based on MSAP data, percentages of methylat- 
ed, hemimethylated, and unmethylated loci were deter- 
mined for each species and for a combination of 
S. microdontum and S. kurtzianum (mid-parental value). 
ANOVAs were performed to compare fragments per indi- 
vidual and methylation percentages between species using 
STATGRAPHICS Centurion XVI software. 

Population structure 

Multivariate analyses were carried out separately for 
genetic and epigenetic data sets. PCoA and cluster analy- 
ses were performed to determine relationships among 
individuals based on genetic and epigenetic polymor- 
phism. Similarity matrices were generated based on the 
Dice coefficient, and the UPGMA linkage method was 
used for dendrogram construction. Also, the Mantel test 
was carried out to contrast genetic, epigenetic, and mor- 
phological matrices. For this, matrices containing shared 
individuals were constructed and similarity matrices, and 
tree topologies were compared using 1000 permutations. 
PCoA, cluster analysis, and the Mantel test were 
performed with NTSYS-pc 2.11, and for bootstrapping 
analysis, the WinBoot program was used (Yap and Nelson 
1996). Genetic and epigenetic distance matrices were 
obtained based on the corresponding Dice similarity 



matrices. For this, each similarity value was transformed 
into distance value (d = 1-s), and a mean distance was 
calculated for each species. ANOVAs were carried out to 
compare genetic and epigenetic variability among species 
using STATGRAPHICS Centurion XVI software. 

Bayesian clustering 

A model-based clustering method was used for inferring 
population structure using allele frequencies. We used 
the program STRUCTURE 2.3 (Pritchard et al. 2000) to 
analyze AFLP and MSAP data sets. Dominant marker 
data were entered by coding both alleles as (1) when 
the fragment was present and both as (0) when the frag- 
ment was absent, specifying (0) as a recessive allele for 
all loci as described by Falush et al. (2007). Estimates for 
the log likelihood were obtained using the admixture 
model and the assumption that the allele frequencies are 
correlated. The log likelihood was estimated from 10 
replicate runs at each K, ranging from K = 1 to K = 6. 
For each run, we used a burn-in of 25,000 cycles and a 
data run of 100,000 cycles. For selecting the optimal 
value of K, similarity coefficients (Rosenberg et al. 2002) 
and AK values (Evanno et al. 2005) were calculated 
using STRUCTURE-SUM (Ehrich et al. 2007) with the R 
software 2.15.1. 

Detection of epiloci under selection 

The MSAP data set was analyzed to identify candidate 
epiloci under natural selection using the Bayesian method 
proposed by Foil and Gaggiotti (2008) and implemented 
in the program BayeScan 2.1 (http://cmpg.unibe.ch/soft- 
ware/bayescan). Assuming that allele frequencies within 
populations follow a Dirichlet distribution, BayeScan 
calculates locus-population F sx coefficients. Selection is 
introduced by decomposing locus population F sx coeffi- 
cients into a population-specific component ft (beta), 
shared by all loci, and a locus-specific component a 
(alpha), shared by all the populations using a logistic 
regression. Departure from neutrality at a given locus is 
assumed when the locus-specific component is necessary 
to explain the observed pattern of diversity (y. signifi- 
cantly different from 0). A positive value of alpha suggests 
diversifying selection, whereas negative values suggest bal- 
ancing or purifying selection. This leads to two alternative 
models for each locus, including or not the alpha compo- 
nent to model selection. BayeScan implements a revers- 
ible-jump MCMC algorithm to estimate the posterior 
probability of each one of these models. Posterior odds 
(PO) are simply the ratio of posterior probabilities and 
indicate how more likely the model with selection is com- 
pared with the neutral model (Foil and Gaggiotti 2008). 
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We used 20 pilot runs of 5000 iterations to adjust 
acceptance rate for each parameter. A burn-in of 50,000 
iterations was employed to attain convergence before 
starting the sampling. The sample size was set to 5000 
with a thinning interval of 10, resulting in a total chain 
length of 100,000 iterations. The loci were ranked accord- 
ing to their estimated posterior probability, and all loci 
with a value greater than 0.91 were retained as outliers. 
This corresponds to log 10 (PO) >1, which provides deci- 
sive support for acceptation of the model according to 
Jeffreys' scale of evidence (Jeffreys 1961). 

Isolation, sequence, and bioinformatic 
analysis of MSAP fragments 

Fragments of interest were excised from the gel and 
heated at 95°C for 20 min in 100 /(L of ultrapure water. 
The eluted DNA fragment was re-amplified with the same 
primers as the MSAP reaction. The PCR protocol was as 
follows: 1 cycle at 94°C for 3 min, 35 cycles at 94°C for 
1 min, 58°C for 30 s, and 72°C for 1 min, and a final 
extension step at 72°C for 20 min. The PCR products 
were resolved in a 1.2% agarose gel, sliced, and purified 
with Wizard® SV Gel and PCR Clean-Up System 
(Promega) following the manufacturer's protocol. Direct 
sequencing of the fragments was performed by Sanger 
technology in an Applied Biosystems Genetic Analyzer 
3130XL. Sequences were edited and analyzed with BioEdit 
7.0.9 (Hall 1999). 

Analysis of DNA similarity was carried out using the 
BLASTn packages implemented at The Institute for Geno- 
mic Research (TIGR; http://www.jcvi.org/; http://blast.jcvi. 
org/euk-blast/plantta_blast.cgi). 

Results 

Collection of plant material 

We collected plant material from the hybridization area, in 
the locality of Guanchin (La Rioja), in two opportunities. 
In the 2008 expedition, we found that the most common 
species was the hybrid S. x rechei. Several populations of 
this species were identified growing in modified environ- 
ments such as backyards, roadsides, and under walnut 
(Juglans sp.) and horticultural orchards. The parental spe- 
cies, S. kurtzianum and S. microdontum, were not present 
in these disturbed areas. We also explored natural environ- 
ments for the presence of the three species. In the gorge of 
the Pismanta River, along a transect approximately 1 km 
long, we sampled six populations of S. x rechei (Table 1). 
No populations of S. kurtzianum were observed, and we 
were able to found only one population of S. microdontum 
(GV07). Then, we explored several gorges in the road to 
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Figure 4. Principal coordinates analysis and minimum spanning tree 
of Solarium kurtzianum, Solarium microdontum, and 
Solanum x rechei based on 30 morphological characters. 

Mina de Oro, and only one population of S. x rechei 
was sampled (GV12). We collected four populations of 
S. kurtzianum separated from the area of hybridization by 
a distance of about of 5 km (CG) and 50 km (LR1, LR2, 
and LR3). In 2010, we expanded the collection area along 
the gorge of the Pismanta River in a transect of about two 
kilometers long. Again, we found that the species with 
greater distribution was the hybrid S. x rechei and sam- 
pled one more population (CMR). In addition, we found 
and sampled one additional population of S. microdontum 
(CMM) and one population of S. kurtzianum (CMK) 
within the gorge of Pismanta River. 

Morphological variation and floral 
abnormalities in 5. x rechei 

We measured 30 qualitative and quantitative characters in 
three plants per accession once flowering had started 
(Table 2) to test for hybridity based on morphological 
intermediacy. There was superimposition in four (14, 19, 
22, and 25) of 30 characters among the three species 
(Table 2). The parental species, S. kurtzianum and 
S. microdontum, showed no significant differences in four 
(11, 18, 20, and 30) additional characters. 

For the character count procedure, only those that 
significantly differed between S. kurtzianum and S. micr- 
odontum were analyzed (22 of 30). Seven characters 
resulted intermediate and 15 resulted non-intermediate 
(Table 2); thus, the one-sided test of intermediate ver- 
sus non-intermediate characters was not significant 
(P > 0.05), and the hypothesis of divergence could not be 
rejected. Six of the non-intermediate characters were 
transgressive in S. x rechei beyond the parental species. 
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Additionally, the association between the three species 
was determined by a multivariate analysis based on the 30 
morphological characters. Figure 4 shows a plot of the 
first two coordinates of the PCoA, which accounted for 
more than 60% of the variability. Individuals identified as 
S. microdontum and S. kurtzianum by their morphotype 
were all gathered into two well-defined groups. On the 
other hand, individuals classified as S. x rechei formed 
three different groups: two located close to each other 
and between S. kurtzianum and S. microdontum, and a 
third group overlapped with S. kurtzianum group. In the 
cluster analysis, S. x rechei individuals were distributed 
in a main cluster and in a second cluster mixed with 
S. kurtzianum individuals, while S. microdontum formed a 
separated group (Fig. Sla). 

Several S. x rechei accessions showed floral abnormali- 
ties (Table 1) that ranged from slight changes in the cor- 
olla form and symmetry to severe alterations such as 
rudimentary stamens and petals; homeotic changes like 
staminoid petals and carpelloid stamens (Fig. S2). 

AFLP and MSAP analyses 

The AFLP analysis with six primer combinations pro- 
duced 832 fragments for the 81 individuals analyzed, 98% 

(A) S. kurtzianum S. microdontum 




S. x rechei 



(546) 

(B) S. kurtzianum S. microdontum 




S. x rechei 



(256) 

Figure 5. Venn diagrams showing the distribution of AFLP (A) and 
MSAP (B) fragments in Solanum kurtzianum, Solarium microdontum, 
and Solanum x rechei. In parenthesis, total number of different 
fragments (loci) per species. 



of which were polymorphic. The error rate based on rep- 
licated samples was 3.2%. The total number of fragments 
obtained was similar for each species, and the mean num- 
ber of fragments per individual showed no significant 
differences among species (Fig. S3a). Of the 276 frag- 
ments shared by the parental species, 33 were absent in 
the hybrid S. x rechei. On the other hand, 47 fragments 
were found in S. x rechei, but not in the parental species. 
The distribution of fragments among the progenitor spe- 
cies and the hybrid is shown in Figure 5A. 

In the MSAP analysis, five primer combinations were 
used and a total of 231 loci were analyzed in 77 individu- 
als. The error rate based on replicated samples was 3.4%. 
Solanum x rechei presented more fragments than the 
other two species, and the number of fragments per indi- 
vidual was significantly higher for S. x rechei than for 
S. kurtzianum (Fig. S3b). Thirty- two loci were scored as 
methylation-insensitive. The other 199 loci (86%) were 
considered methylation-sensitive and used to determine 
the epigenetic diversity. The percentage of unmethylated 
loci for S. x rechei (Mean = 48.0%) was significantly 
(P < 0.01) higher than that for the mid-parental value 
(Mean = 44.7%). The percentage of methylated loci for 
S. x rechei (Mean = 40.5%) was significantly (P < 0.05) 
lower than that for the mid-parental value 
(Mean = 43.9%). The difference in the hemi-methylation 
level was not statistically significant (P > 0.05; Fig. 6). 

The Venn diagram in Figure 5B shows methylation 
patterns common to two or all three species and unique 
patterns for each species. The hybrid S. x rechei pre- 
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Figure 6. Percentage of methylated, hemimethylated, and 
unmethylated loci based on MSAP data for Solanum kurtzianum (ktz), 
Solanum microdontum (mcd), Solanum x rechei (rch), and a mid- 
parental value. Significance levels are **P < 0.01 , *P < 0.05. 
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Figure 7. Genetic and epigenetic distances (Dice coefficient) of 
Solanum kurtzianum (ktz), Solanum microdontum (mcd), and 
Solanum x rechei (rch) based on AFLP and MSAP markers, 
respectively. Different letters indicate statistically significant difference 
at P< 0.01. 

sented 43 exclusive patterns, not observed in the parental 
species, representing 12% of all patterns detected. Among 
them, 29 (8%) are unequivocally methylation changes; as 
for the same loci, the parental species showed the pres- 
ence of the fragments, but with a different pattern respect 
to the hybrid. For the other 14 loci (4%), the parental 
species showed no fragments, which could indicate both a 
change in the methylation status in the hybrid (if the sites 
where hypermethylated in the parents) or a genetic 
variation. 

Genetic and epigenetic structure 

Genetic and epigenetic distances among individuals of the 
same species were used to calculate the genetic and epige- 
netic variability of each species. No significant differences 
in variability (genetic or epigenetic) were found between 
species (P > 0.05), although the comparison of genetic 
against epigenetic distances within species was statistically 
significant for the three species (P < 0.01; Fig. 7). 

The two-first principal coordinates of the PCoA with 
the AFLP and the MSAP data sets (Fig. 8A,B) accounted 
for 32.8% and 24% of the variance, respectively. Geneti- 
cally, S. kurtzianum group showed a greater aggregation 
of the individuals than epigenetically. In the first case, all 
individuals were clustered together, while in the second 
case, genotypes from La Rioja separated from those 
accessions from Mendoza. Solanum microdontum behaved 
similarly in both studies, where plants collected in La 
Rioja province and the accessions from Catamarca and Sal- 
ta provinces provided by the PFGB established two distinct 
subgroups. Solanum x rechei individuals in both plots 
were located between the parental species; although in the 
MSAP plot, they were considerably closer to S. kurtzia- 
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Figure 8. Principal coordinates analyses and minimum spanning tree 
of Solanum kurtzianum, Solanum microdontum, and 
Solanum x rechei individuals based on nuclear AFLP markers (A) and 
MSAP patterns (B). 



Table 4. Pairwise comparison of MSAP, AFLP, and morphology matri- 
ces with the Mantel test. 





MSAP 


AFLP 


Morphology 


MSAP 
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-0.628 


AFLP 
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-0.5S0 (56) 
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Below the diagonal: comparison between similarity matrices. Above 
the diagonal: comparison of tree topologies transformed to matrix val- 
ues through cophenetic correlation. Numbers in parentheses indicate 
the number of accessions common to each comparison. 0 = no corre- 
spondence, 1 = perfect correspondence. P < 0.01 was obtained for 
all comparisons based on 1000 permutations test. 

num. Additionally, in the AFLP analysis, the hybrid was 
dispersed in many subgroups, which corresponded to the 
different sampled natural populations, while in the MSAP 
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Figure 9. Bar plots showing genetic and epigenetic structure of Solanum kurtzianum, Solanum microdontum, and Solanum x rechei populations 
based on AFLP and MSAP data, respectively. Each accession is represented by a vertical line partitioned in K colors that represent the proportion 
of membership to each cluster. For AFLP, K = 2 and 3 and for MSAP, K = 3 are shown. 



analysis, the group was more condensed except for popula- 
tions GV02 and CMR. In the cluster analysis, individuals 
were grouped similarly as in the PCoA (Fig. Slb,c). 

The Mantel tests (Table 4) showed a higher correlation 
of morphological data with MSAP (r = 0.55, P < 0.01) 
than with AFLP (r = 0.38, P < 0.01). MSAP matrix holds 
a greater correlation with AFLP (r = 0.78) given that both 
are methodologically related (i.e., for two individuals to 
have the same methylation state at any given locus, they 
must first share the locus). 

Population structure 

To choose the optimal number of groups for the STRUC- 
TURE analysis, 10 runs were performed for each value of 
K from 1 to 6. For the AFLP data set, K = 2 and 3 pre- 
sented similarity coefficients close to 1, which means that 
individuals were repeatedly assigned to the same groups 
in all runs. This behavior was also observed in the bar 
plots. With MSAP data, the similarity coefficient was 1 
only for K = 3. Furthermore, the AK values suggested 
fT=2 as the optimal for AFLP and K = 3 for MSAP 
(Fig. S4). 

In Figure 9, the populations genetic structure for K = 2 
and 3 and epigenetic structure for K = 3 are represented. 
Each vertical line corresponds to an individual, and each 
color symbolizes the relative percentage of membership to 
each cluster. The parental species, S. kurtzianum and 
S. microdontum, corresponded to one cluster each 
and showed negligible admixture levels for both AFLP 
and MSAP data. In contrast, S. x rechei differed greatly 
between its genetic and epigenetic structure. For AFLP 
data, when K = 2 was tested, and all S. x rechei individu- 
als presented certain percentage of membership to both 
parental clusters, and for K = 3, some individuals were 



assigned to a third group, while the rest presented high 
levels of admixture with the parental groups. For MSAP 
data, most individuals of S. x rechei belonged to their 
own cluster. 

Detection of epiloci under selection 

BayeScan analysis identified 24 MSAP epiloci under selec- 
tion (or linked to selected loci), with "strong evidence" 
(log 10 (PO) > 1) on the Jeffreys' scale. Fifteen of these 
markers were pointed out "very strongly" to be under 
selection (log 10 (PO) > 1.5), and for seven of them, evi- 
dence was considered "decisive" (log 10 (PO) > 2; Fig. S5). 
All these epiloci were characterized by high F ST coeffi- 
cients and positive i values and are therefore interpreted 
as having been affected by divergent selection. 

Four MSAP fragments with the highest F ST coefficients 
were isolated and sequenced. Two sequences were of high 
quality and were used for BLAST search on the TIGR 
database. One MSAP fragment of 144 bp was similar 
(E value 1.6 x 10~ 13 ) to an EST (Expressed Sequence 
Tag) from Solanum chacoense (GeneBank Accession 
DN978520). The other fragment of 125 bp was similar (E 
value 4.9 x 10~ 7 ) to a sequence of Solanum lycopersicum 
(GeneBank Accession TA55928_4081) annotated as a 
pararetrovirus-like sequence reported in tobacco. 

Discussion 

In the present work, we analyzed the genetic and epige- 
netic variation generated in the recently formed natural 
potato hybrid S. x rechei based on a genome-wide 
screening method. The hybrid nature of S. x rechei had 
been analyzed using morphological and molecular 
(RFLPs) approaches by Okada and Hawkes (1978) and 
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Clausen and Spooner (1998), respectively. However, 
Clausen and Spooner (1998) used for the RFLP analysis 
the diploid accessions of the parental species S. microdon- 
tum from non sympatric regions northern of the hybrid- 
ization area in the province of La Rioja, Argentina. Here, 
we studied plants of S. microdontum from two diploid 
accessions (accession number 4631 and 5902 from Cata- 
marca and Salta provinces, respectively) used by Clausen 
and Spooner (1998), and in addition, we collected trip- 
loid accessions (M. S. Ferrer, C. F. Marfil, and R. W. 
Masuelli, unpublished results) of S. microdontum from 
Guanchm (La Rioja province) growing sympatrically with 
S. x rechei. In S. microdontum, PcoA analyses of genetic 
data showed intraspecific discontinuities, separating 
populations from La Rioja province in one group and 
the populations from Catamarca and Salta provinces in 
another group. However, these differences were not sup- 
ported in the Bayesian STRUCTURE analysis. This result 
confirms that S. microdontum material used by Clausen 
and Spooner (1998) was representative of the parental 
populations. 

As described by Marfil et al. (2009) for a different 
potato hybrid {Solanum ruiz-lealii), S. x rechei presented 
genotypes with different floral abnormalities, suggesting 
that recently formed hybrids may suffer genetic and 
epigenetic instability as a result of the interspecific 
hybridization event. 

Besides the general additivity of the hybrid respect to 
the parental genomes, gain (5.6%) and loss (4%) of AFLP 
fragments were observed in the hybrid, possibly as a 
result of genome rearrangements induced by the interspe- 
cific hybridization. Similar results were described in inter- 
specific hybrids of wheat (Liu et al. 1998), Spartina spp. 
(Salmon et al. 2005) and Solanum (Marfil et al. 2006). In 
the hybrid, the epigenetic modification was more signifi- 
cant than the genetic changes, showing twice (12%) the 
number of novel fragments in MSAP than in the AFLP 
analysis. The MSAP loci analyzed showed a general deme- 
thylation compared with the parental genomes. Plants of 
S. kurtzianum and S. microdontum had 43% and 45% 
methylated loci, respectively, while 40% of the loci from 
S. x rechei were methylated. Zhao et al. (2007), com- 
pared relative levels of cytosine methylation at the CCGG 
sites and observed demethylation in maize hybrids respect 
to their parental inbred lines. Previous studies have 
showed that chemical demethylation of the potato hybrid 
S. ruiz-lealii induced floral abnormalities similar to the 
phenotypes observed in natural populations (Marfil et al. 
2009). In addition, early flowering, changes in leaf 
morphology, and up-regulation of miRNAs were observed 
in plants treated with the demethylating agent AzaC 
(Marfil et al. 2012). Additional evidence of genomewide 
demethylation associated with runaway transposition of 



retroelements comes from interspecific hybridization 
between two wallaby species (Waugh O'Neill et al. 1998) 
and from newly synthetic Arabidopsis allopolyploids that 
showed loss of cytosine methylation in retrotransposon- 
related elements (Comai et al. 2000). 

Our study showed a high level of genetic and epigenetic 
polymorphism for the three taxa analyzed, where 98% 
and 99% of the AFLP and MSAP loci were polymorphic, 
respectively. For each of the studied taxa, the epigenetic 
variability was significantly higher than the genetic vari- 
ability, based on Dice coefficient (Fig. 7). PcoA analyses 
of morphological, genetic, and epigenetic information 
differentiated the three taxa and grouped the hybrids in 
intermediate position between the two parental species 
(Figs. 4, 8). The Bayesian STRUCTURE analysis of the 
genetic data showed that S. x rechei accessions, except 
for accession GV08 and one individual from CMR and 
GV11, were assigned to a cluster with posterior probabil- 
ity lower than 0.5 (Fig. 9). These results reveal that 
S. x rechei individuals have different levels of genetic 
admixture, indicating a recent hybridization and intro- 
gression between S. kurtzianum and S. microdontum. 
Interestingly, pairwise comparison analysis among the 
three matrices showed that the morphological data corre- 
lated better (r = 0.550, P < 0.01) with the epigenetic than 
the genetic information (r = 0.383, P < 0.01; Table 4). 
These results suggest that, in a recently formed hybrid, 
the epigenetic variation influences the phenotypic out- 
come more than the genetic one. This conclusion was 
supported by the Bayesian STRUCTURE analysis of the 
genetic and epigenetic data. The genetic data showed the 
hybrid nature of S. x rechei with an estimated number of 
groups of 2 {K = 2) with individuals displaying admixture 
between S. kurtzianum and S. microdontum. Instead, the 
epigenetic information clearly assigned S. x rechei indi- 
viduals to a third cluster (K = 3) with high probability 
(Fig. 9 and Fig. S4), with seven individuals of 30 with 
more than 10% of admixture. These results indicate that 
an epigenetic pattern was rapidly established in the 
hybrids differentiating them from the parental taxa. In 
allopolyploid orchids, Dactylorhiza, Paun et al. (2010) 
showed that the methylation patterns explain better the 
gene expression differences between allopolyploids than 
the genetic information. It is noteworthy to show that in 
S. x rechei the hybridization event rapidly established a 
novel epigenetic pattern, different from the parental spe- 
cies, which could modulate the regulatory networks and 
persist stably over many generations possibly through 
clonal propagation by tubers. According to our results, 
the adaptation of S. x rechei to grow as a weed in dis- 
turbed environments, such as along roadsides or under 
trees in orchards (Okada and Hawkes 1978), is better 
explained by the rapid setup, after the hybridization 



3776 



© 201 3 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



N. Cara ef al. 



Epigenetics and Hybridization in Solarium 



event, of an epigenetic pattern influencing phenotypic 
plasticity and adaptation to environmental changes. Baye- 
Scan analysis identified 24 epialleles under selection, seven 
of them with high F sx values and therefore interpreted as 
having been affected probably by divergent selection. Nev- 
ertheless, it is possible that these epialleles may have 
arisen as a consequence of the hybridization event and 
have been maintained by clonal propagation. One of the 
sequenced epialleles was related to an endogenous para- 
retrovirus sequence reported in tobacco, possibly confer- 
ring heritable virus resistance (Mette et al. 2002). The 
authors found that the pararetrovirus sequence was per- 
sistently methylated in CNG and CG sites, and the meth- 
ylation patterns were preserved in two tobacco species, 
the allotetraploid Nicotiana tabacum and the diploid 
N. sylvestris. These results suggest that the epiallele was 
selected in both species and conferred advantageous traits 
to tobacco and N. sylvestris. It is likely that the same 
epiallele was under selection in the potato species studied. 

Changes in chromosome number could play a role in 
the origin of the variability naturally found in the hybrid 
populations. As the hybrids were probably originated by 
the hybridization of the diploid S. kurtzianum and trip- 
loid cytotypes of S. microdontum, some hybrid genotypes 
are likely to be aneuploids. However, Okada and Hawkes 
(1978) found only euploid (both diploid and triploids) 
hybrid cytotypes in nature. We are currently analyzing in 
natural and synthetic hybrid plants the occurrence of 
aneuploids and the effect on plant phenotypes. Prelimin- 
ary results indicate that some hybrid plants are aneup- 
loids (M. S. Ferrer, C. F. Marfil and R. W. Masuelli, 
unpublished results). 

Concluding Remarks 

Hybridization is an important event in the evolution of 
the potato species, and several authors suggested that 
more than 10% of the Solarium species from section 
Petota are the result of hybridization. There is also an 
important amount of phenotypic plasticity in this section 
(Correll 1962; Hawkes and Hjerting 1969). Our results 
show that after hybridization, a novel epigenetic pattern 
is rapidly established and maintained through several 
generations allowing the adaptation of hybrids to new 
environments. The importance of epialleles in the evolu- 
tionary process is intriguing. We found that for the three 
species, the methylation polymorphism detected was 
higher than nuclear DNA variability; and in addition, a 
better correlation between epigenetic patterns and pheno- 
typic characters was observed. According to Kalisz and 
Purugganan (2004), given that epiallelic phenotypes are 
more flexible than alleles produced by genetic mutation, 
it is possible that phenotypes produced by epialleles 



could explore novel environments and drive the adapta- 
tion of other genes in the short term. We strongly 
suggest that epigenetic variation, together with genetic 
variation, plays a relevant role in the evolution of the 
Solarium species section Petota and that the analysis of 
the epigenetic variation may shed some light on the intri- 
cate taxonomy and systematic relationships among 
species of this section. 

Data Accessibility 
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Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Figure SI. Cluster analyses of Solatium kurtzianum, Sola- 
tium tnicrodontum and Solatium x rechei based on mor- 
phological (a), AFLP (b) and MSAP (c) data. For the 
AFLP and MSAP dendrograms bootstrap values are 
shown above the branches. 

Figure S2. Flower phenotypes of Solatium x rechei indi- 
viduals. 

Figure S3. Total number of fragments per species and 
mean number of fragments per individual obtained by 
AFLP (a) and MSAP (b) techniques. 
Figure S4. Similarity coefficients and Delta K values for 
AFLP and MSAP datasets. 

Figure S5. Scan to identify epiloci subject to selection 
using BayeScan. 
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